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Abstract 

Wounds in the embryo show a remarkable ability to heal quickly without leaving a scar. Previous 
studies have found that an actomyosin purse string forms around the perimeter of the wound to help 
close it over the course of several dozens of minutes. Here, we report on an even faster mechanism 
which closes an embryonic wound by more than 50% within the first 30 seconds. Circular and 
elliptical wounds (~100 fj,m in size) were made in the blastoderm of early chick embryos and 
allowed to heal, with wound area and shape characterized as functions of time. The closure rate 
displayed a biphasic behavior, with rapid constriction lasting about a minute, followed by a period 
of more gradual closure to complete healing. Fluorescent staining revealed that, within seconds 
after wounding, actin and myosin-II increased in intensity at cell borders within a relatively thick 
ring (several cells wide) around the wound. By one minute post wounding, this structure gave 
way to a narrow actomyosin band at the wound border, consistent with a "purse string" healing 
mechanism. We hypothesized that contraction of the thick ring is responsible for the initial, rapid 
phase of wound closure, while formation and contraction of a thinner purse string drives the slow 
phase. To test this hypothesis, we implemented both contractile mechanisms in a finite-element 
model. Comparing model and experimental results for circular and elliptical wounds supports our 
hypothesis and indicates that the initial cellular contraction is relatively isotropic in nature. The 
model also explains an otherwise counterintuitive observation, that some elliptical wounds open 
briefly before going on to heal. The results suggest that rapid cellular contraction works together 
with slower contraction of an actomyosin ring to heal embryonic wounds. 

Keywords: biomechanics, chick embryo, epithelium, finite elements, epithelial morphogenesis, 
mechanobiology, growth 



1. Introduction 



Embryonic epithelia display a remarkable ability to repair wounds, and use a series of redundant 
mechanisms to do so. Shortly after a wound is made, an actomyosin cable forms at the wound edge 



and contracts, drawing the wound closed by a purse string mechanism (Jacinto et al. , 2001 ; Martin 



and Lewis, 1992; Redd et al. 



2004). Filopodia and lamellipodia also frequently form to draw in 



apposing edges like a zipper (Jacinto et al. 2000 Wood et al. 2002). Both these mechanisms 



require the polymerization of new actin structures and can form relatively quickly, over a period of 
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minutes. Here, we report on a cell contraction mechanism that precedes these other mechanisms 
and operates on an even faster time scale. Together, these mechanisms constitute three phases for 
embryonic wound healing. 



Much of the cellular machinery used in wound healing also is used during morphogenesis (Martin 



and Parkhurst, 2004 Sonnemann and Bement, 2011 Wood et al. 2002). In Drosophila dorsal 



closure, for example, two epithelial sheets are drawn together and fuse in a process remarkably 
similar to the healing of embryonic wounds, with both the purse string and zippering mechanisms 
playing a significant role (Hutson et al. 2003). The mechanisms of embryonic wound healing 



are significantly different from adult wounds, which tend to involve crawling fibroblasts and an 
inflammatory response that can lead to scarring. Embryonic wounds, by contrast, heal relatively 
quickly and do not form scars, and thus may offer clinical insight to improving adult wound 



healing (Redd et al. 2004; Sonnemann and Bement 2011). 



Here, we studied the mechanics of wound healing in the blastoderm of early chick embryos. 
We found that a rapid healing phase partly closes the wound within tens of seconds, followed by 
a slower healing phase which continues over the course of minutes. Using fluorescence microscopy 
and finite-element modeling to characterize the healing mechanisms, we found that contraction 
within a relatively thick ring of cells around the wound is responsible for the initial phase of wound 
closure, while assembly and contraction of an actomyosin cable ("purse string") accounts for the 
second, slower phase. A third phase, not studied here, is "zippering" as filopodia pull wound edges 
together to fuse when they approach each other. 



2. Background 



The chick embryo has been a popular model for studies of wound healing (Bortier et al. 1993 



1977; 



Martin and Lewis 



1992: 



Brock et al.[ |1996[ |England and Cowper[ |1977| |Lawson and England} |1998[ |Mareel and Vakaet 



Stanisstreet et al. 



1980). The early embryo consists of a planar 



blastoderm, two to three cell layers thick, which is held under tension and attached to a substrate 



(the vitelline membrane) only at its periphery (Bellairs et al. 1967 Bortier et al. 1993 New 



1959). All forces required to close a wound are generated within the blastoderm itself. 

A critical initial signal of a wounding event is the influx of extracellular calcium (Ca 2+ ) into 
the perforated cells, as well as its release from intracellular stores (Benink and Bement[ 2005 



Cordeiro and Jacinto 2013 Woolley and Martin 2000 Xu and Chisholm, 2011). Among other 
roles, Ca 2+ triggers a signaling cascade involving small GTPase molecules which play a central and 
necessary role in the regulation of wound healing, including recruitment of actin and myosin to 



the wound edge, actomyosin assembly and contraction, and the formation of filopodia (Bement 



eFah] [2006] |Brock et alT| [T996] |Clark et al.] [20091 |Sonnemann and Bemehtj [20TT| |Wood et aL 



2002). A number of different wound healing mechanisms occur in various model organisms. Ac- 



tomyosin purse strings are widely observed in both single-cell and multicellular wounds (Clark 



Kiehart, 1999; Martin and Lewis, 1992). Epithelial cells crawling over a mesenchymal 
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the wounds were cut completely through the blastoderm, and there is no substrate to aid healing. 

Published work on the healing of elliptical-shaped embryonic wounds is inconclusive on how 
the shape of the wound changes as it heals. Elongated wounds in Drosophila embryos maintain 



an approximately constant aspect ratio (AR) as they close (Hutson et al. 2003). Small elliptical 



wounds in Xenopus oocytes (frog eggs) round as they heal (Mandato and Bement, 2001). Larger 
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superficial rectangular wounds in multicellular Xenopus embryos with an underlying mesenchyme, 



however, become more elongated during healing (Davidson et al. 2002]). 

In contrast to wound healing in adult tissue, which operates by markedly different processes, 



modeling of embryonic wound healing has received relatively little attention ( Murray , 2003 , Olsen 



et al. , 1995 Wyczalkowski et al. 2012 ) . Much of the published work is based on the mechanochem- 



Sherratt et al. 



ical model of Murray and Oster (1984) for the morphogenesis of epithelial sheets 
( 1992|) (see also Murray||2003[) considered the quasi-static response of such a sheet to wounding 
and Sadovsky and Wan| ( |2007| ) extended this model to incorporate wound closure. Previously we 
presented a model for a circular wound in which the healing response is governed by a stretch- 



activated morphomechanical feedback law (Taber 2009). Nagai and Honda (2009) developed a 



vertex dynamics model where minimization of interfacial energy leads to the closing of a wound. 



Finally, Hutson et al. (2003) (see also Layton et al. 2009) developed a model incorporating an 



actomyosin ring and filopodial zippering for the related problem of dorsal closure in Drosophila. 

Tissue-scale deformations are driven by cell division, cell intercalation, cytoskeletal contraction, 
or some combination of these and other processes. The present analysis is based on our tissue- level 



theory for epithelial morphogenesis ( Taber , 2009 ) , in which these processes are simulated by active 



changes in the local zero-stress configuration. Here, we consider only active contraction of tissue, 
as simulated by negative growth, and do not consider the contributions of filopodia. 



3. Experimental Methods 

3.1. Embryo Preparation and Culture 

Fertilized white Leghorn chicken eggs were incubated at 37°C for 20-22 hours (of a 21-day 
incubation period) in a humidified, forced draft incubator to yield embryos at approximately 



Hamburger-Hamilton (HH) stage 4 (Hamburger and Hamilton 1951). Whole embryos were har- 



vested from the eggs using a filter paper carrier method (Voronov and Taber 2002) which preserves 
the stresses normally present in the tissue. 

At this early stage of development, the embryonic blastoderm is organized as a nearly flat 
sheet in a state of approximately uniform, isotropic tension (Varner et al. 2010) (Fig. la). The 



blastoderm is connected to the vitelline membrane only at its outer periphery, and consists of two 
epithelial cell layers, the endoderm and ectoderm. A third cell layer, the mesoderm, consists of 
migratory mesenchymal cells and spreads between the endoderm and ectoderm, although at this 
stage of development its extent is limited (Bortier et al. 1993 England and Wakely, 1977). While 



the endoderm and ectoderm have important differences, wounds passing through both layers heal 
in the same manner as ectoderm-only wounds ( Stanisstreet et al. 1980), and in characterizing 



wound healing dynamics we do not distinguish between the layers. 

After harvesting, embryos were placed atop a 2% agar/98% phosphate buffered saline (PBS) 
gel in Delta T Dishes (Bioptechs, Butler, PA), submerged under a thin layer of liquid culture 
media ( Voronov and Taber[ 2002), and incubated for at least one hour. Prior to wounding, the 
embryos were removed from the incubation chamber and allowed to equilibrate in the Delta T4 
Culture Dish Controller under the same incubation conditions. 



3.2. Wounding and Imaging 

Circular and linear incision wounds were cut completely through the blastoderm, 
wounds used for dynamic imaging were positioned as illustrated in Fig. la 



Circular 
Elliptical wounds, 
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Figure 1: Early stage (HH4) chick embryos with circular and elliptical wounds, ventral side up and anterior 
to the left. Hensen's node marked by "H" . Scale bars 500/im. |(a)| Circular wounds (numbered 1-6) several 
seconds after being created with a circular punch (inset). |(b)| Two elliptical wounds (indicated by *) created 
with microscalpel (inset). 



as well as circular wounds used for fluorescence imaging, were positioned anterior and lateral to 



Hensen's node (as in Fig. lb), where the mesoderm is absent (England and Wakely, 1977). The 



location of the wounds generally did not affect their behavior. 

To create circular wounds, a microforged glass micropipette (42 //m inner, 57 ^m outer radii) 
held in a micromanipulator was used to punch through the blastoderm and excise circular plugs of 



tissue of uniform size (Fig. la)(Varner and Taber, 2010). Linear incisions approximately 200 fxm in 



length were made by drawing a handheld microscalpel (Conrad et al. 1993) across the blastoderm 
surface. Tension in the membrane caused these wounds to open into an initial elliptical shape 



(Fig. lb). In both cases, wounds were concealed by the wounding apparatus while they were made, 
with the first unobscured frame visible 2-4 seconds after wounding. 

Wounds were created and allowed to heal under incubation conditions, and were observed for 
the entire duration of the experiment by a combination of video and still digital imaging. Wound 
area was determined by thresholding each image to isolate the wound, with manual correction as 



necessary, using ImageJ (Abramoff et al. 2004). The ratio of major to minor radius of elliptical 



wounds (the wound aspect ratio, AR) was obtained by fitting the isolated wound shape with an 
ellipse of equivalent area and second moment of inertia. 

Fluorescence staining was used to visualize the actin cytoskeleton and phosophorylated myosin 



II light chain (pMLC), as described previously (Filas et al. 2012). Wounded embryos were fixed 
at a given time after wounding by injecting 37% formaldehyde into the media at a 1:10 volume 
ratio. F-actin was visualized with rhodamine phalloidin, while immunofluorescence staining was 
used for pMLC. Following staining, the endoderm and ectoderm in the vicinity of the wound were 
separated using drawn glass capillary tube needles. Small pieces of both germ layers were then 
excised and mounted for fluorescence imaging. 



4. Experimental Results 

We analyzed 52 circular wounds from 15 embryos and 13 elliptical wounds from 10 embryos. 
Only wounds which closed completely, or nearly so, during the period of observation were included 
in this analysis. Nearly all (62/65) wounds closed relatively quickly and typically healed fully within 
ten minutes to two hours, depending on wound size as well as wound location and stage of embryonic 
development. Three circular wounds initially closed partly, then expanded for several minutes 
before going on to heal, and a minority of wounds never healed and were not analyzed. Relatively 



wide variations in healing time for such wounds have been previously reported ( Stanisstreet et al. 
19801). 



Our analysis focuses on the first ten minutes after wounding. Most of the healing, as measured 
by wound area, took place during this time, and optical characterization of wound area and shape 
is most reliable during this period, as discussed below. Subsequent modeling focuses on this time 
period as well. 

4-1. Circular Wound Area 

The experimental wound area is shown as a function of time after wounding for six circular 
wounds (Fig. [2] and Supplemental Video 1). Immediately after formation, wound area decreased 
rapidly for about 30 seconds. After this period the wounds continued to close, but significantly more 
slowly. To minimize variability from differing experimental conditions and embryos, we consider 



here only wounds from a single embryo (shown in Fig. la). However, the qualitative behavior 
shown in Fig. [2] is representative of that observed in nearly all embryos we analyzed; in particular, 
wounds broadly displayed distinct fast and slow phases of healing. 




Figure 2: Experimental wound area vs. time for circular wounds. Experimental data (solid curves) are for 
six wounds from a single representative embryo (pictured in Fi g. [Ta| . Dashed curves indicate best fit of 
Eq. ([I]) for each wound (parameter values in Table |Sl| ) . Panels [a[ |b| show the first 10 minutes and 1 minute 
of wound healing, respectively. Good fit to data obscures dashed curves in panel a. 



We approximate the area versus time relationship for each of the wounds in Fig. [2] using a 
double exponential function of the form 

a(t) = c + ci exp(-t/ri) + c 2 exp(-i/r 2 ), (1) 

which contains two characteristic time constants, n and t 2 . Choosing t\ < r 2 , we identify c\ and 
C2 as the magnitudes of the fast and slow phases of wound healing, respectively. For each wound 
we used nonlinear least squares fitting to obtain values for each of the five parameters in Eq. Q 
(listed in Table SI). The dashed curves in Fig. [2] show the resulting fit, illustrating that the double 



exponential function reproduces trends in the data quite well. 

For all six wounds the fast phase is responsible for the bulk of the decrease in wound area 
(ci ~ 2C2), and the associated time constant is t\ = 12 A ± 2.5 sec. The time constant of the 
slow phase is more variable, but is generally ten or more times greater than n (r 2 = 140 ± 40 
sec). We used the mean parameters and Cj (Table [Slj) to construct a mean wound response. 
Not all wounds closed within ten minutes, an aspect which is captured by the Co term, and which 
accounts for the large variation in c$ . Evaluating Eq. ([I]) for t = with the mean wound parameters 
yields the reference wound area ao = 9130 ^m 2 , corresponding to a radius r w = 54 //m. This initial 
wound size is larger than the pipette inner radius (42 //m), which is the cutting edge. This response 



represents recoil of the wound edge due to blastoderm tension (Varner and Taber 2010) 



The relatively large difference between the time constants T\ and T2 suggests that there are (at 
least) two distinct physical mechanisms operating to close the wound during the first ten minutes. 
The initial phase, associated with time constant t±, results in the bulk of wound closing within ~20- 
30 seconds of wounding (area at 30 seconds is 35% ao). This is followed by a second phase, associated 
with T2, which operates more slowly to close the wound the rest of the way. We hypothesize that 
the second healing phase corresponds to actomyosin cable assembly and contraction, with a time 



scale consistent with published work on the dynamics of actin cable formation ( Abreu-Blanco et al. 

As discussed later, however, the rate of the initial 



2011 


Brock et al. 


1996 


Clark et al. 


2009) 



mechanism appears too rapid for de novo actomyosin assembly. 

4-2. Elliptical Wound Area and Shape 

Following linear cuts, wounds opened into elliptial shapes as quickly as could be observed. The 
area and AR of five representative elliptical wounds, all from different embryos, are plotted in 
Figure [3j The initial areas of the elliptical wounds are more variable than the circular ones, but 
are adequate for a qualitative description of wound healing. 
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Figure 3: Experimental wound area and aspect ratio (AR) vs. time for elliptical wounds. Each wound is 
from a different embryo. ([a[ [b]) Area and AR, respectively, for ten minutes post wounding. (JcJ [d]) The same 
quantities for the first minute. Sampling rate changes correspond to transition from video to still image 



acquisition (see Sec. 3.2 1 



Like circular lesions, elliptical wounds closed at a rapid initial rate, then more gradually until 
fully healed. For some wounds the area initially increased for several seconds before decreasing, 



although this was not observed in all elliptical wounds (Figs. 3a[c). No circular wounds displayed 
such transient gaping behavior. 

Figures 3b|d show the AR of healing elliptical wounds. The AR, obtained from the quotient 
of two fluctuating quantities, is subject to large fluctuations, particularly as wounds become small 



(see also Section 7.2.3). Nevertheless, it is clear the AR tends to stay relatively constant, indicating 
that the wounds neither round up (AR— > 1) nor become slitdike (AR — > 00) as they heal. Notably, 
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in cases where the area of the wound transiently increases immediately after wounding, there is a 
corresponding initial decrease in the AR. 

We reason that, since the same chemo- mechanical processes are involved for both circular and 
elliptical wound healing, differences in their behavior (e.g., initial increase in area) are likely caused 
by geometric effects. 



4-3. Actin and Myosin Staining 

To obtain clues concerning the origin of the forces that close the wounds, we investigated the 
distributions of F-actin and activated myosin (pMLC) at various times after wounding. To visu- 
alize actin distribution more clearly, we separated the two epithelial monolayers, endoderm and 
ectoderm, and imaged them separately (Fig. [4]). These layers display different cell morphologies, 
with the thicker ectoderm consisting of columnar epithelial cells and the thinner endoderm of squa- 



mous cells (Bellairs and Osmond, 2005 England and Lawson, 1993). Nevertheless, the conclusions 



drawn here regarding their healing response are broadly similar. 
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Figure 4: Fluorescence images of circular wounds stained for actin at given times after wounding, 
ectoderm; (|d]jf|) endoderm. Scale bars lOO^m. Panels [b] and [ej [c] and [f] are same wound. Relatively large 
size of wound at later times is an artifact of fixation. 

Tens of seconds after wounding, we observed a broad, diffuse ring of enhanced actin fluorescence 



near the wound in both the ectoderm (Fig. 4a) and the endoderm (Fig. 4d). This intensified 
fluorescence, an indicator of increased F-actin localization, coincided mainly with cell borders. A 
minute after wounding this ring had largely dissipated, except at the wound border where the 
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stain became more intense (Figs. 4b|e ) . By ten minutes actin had formed a continuous, distinct 
supracellular structure along the wound border, characteristic of an actomyosin cable (Figs. 4c|f ). 
For convenience, we denote the initial and subsequent actin regions as the "thick ring" and "thin 
ring," respectively. In addition to the cable structure, finger-like filopodia and leaf-like lamellipodia 



were seen in the endoderm (Figs. 4e, SI), consistent with previous reports (see Section 7.2.3). 

While the boundary of the thick ring was not sharply defined (Fig. [4a| , it was roughly 40 /xm 
in width and of comparable size around both circular and elliptical wounds (data not shown). The 
width of the thin ring was somewhat variable but generally 1-3 /im wide, consistent with published 



reports (Martin and Lewis, 1992). 

To better establish the correlation between actin fluorescence intensity and mechanical force 
production, we also stained circular wounds for phosphorylated myosin II light chain (pMLC), 
which is the activated conformation of the myosin molecule (Lecuit et al. 2011). The results 



(Fig. S2 ) show essentially the same features as actin staining, with a thick ring around the wound 
immediately after wounding and a sharp, condensed cable structure ten minutes later. 

These results suggest that the initial phase of wound healing, associated with the time constant 
n, is driven by rapid cellular contraction within a thick ring of cells around the wound. The second 
phase, associated with the time constant T2, is in turn driven by slower contraction of a relatively 
thin actomyosin cable. 



5. Computational Model for Wound Healing 

To investigate the plausibility of our proposed healing mechanisms, we constructed a finite- 
element model that incorporates active tissue contraction. The blastoderm is modeled as an initially 
homogeneous single-layered membrane in plane stress under isotropic tension <tb (Fig. 5a). The 



material properties associated with this epithelial membrane capture the collective mechanical 
contributions of both the endoderm and ectoderm, including cellular cytoplasm and cortical actin. 
From our experimental results, we propose that wounding triggers two distinct mechanisms: (1) 
initial contraction of cells in a thick ring surrounding the wound, and (2) subsequent assembly and 
contraction of actomyosin fibers in a thin ring at the wound border. 

To incorporate these mechanisms into our model, the membrane is treated as a constrained 
pseudo-elastic mixture consisting of two contractile components: "cells" and "fibers" (Humphrey 



and Rajagopal 2002). Immediately after wounding, cells contract within the thick ring (red area, 
Fig. [5]), while contractile fibers form more gradually in the thin ring, replacing cells in the process 
(blue area, Fig. [5j). The fibers are assumed to be stiffer than cells, created in a prestretched 
configuration, and to contract along their length. Working together, cell contraction, assembly 
of prestretched fibers, and subsequent fiber contraction generate the mechanical forces that close 



the wound. Both circular and elliptical wounds are considered in our model (see Appendix A for 
geometric considerations). 



et al. (1994) for finite volumetric growth ( Ramasubramanian et al. 2006 Ramasubramanian and 



the theory of 


Rodriguez 


Ramasubramanian and 



Taber, 2008), and we consider two models for fiber formation. First, we discuss the "multiple fiber" 
model, which is based on the constrained mixture and evolving natural configurations theory of 
Humphrey and Rajagopal (2002). Next, we describe a simpler "single fiber" model, which is a 



special case of the multiple fiber model and offers numerical efficiency. 
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Figure 5: Schematic of model. |(a)| Geometry of circular wound model in reference configuration bo- The 
blastoderm is modeled as a square single-layered epithelial membrane under isotropic tension erg. |(b)| Detail 
of wound region for circular and elliptical wounds. The thick (cell) ring is shown in red, and the thin (fiber) 
ring in blue. The layers indicate the geometry in the undeformed (-Bo) and reference (bo) configurations, 
illustrating differences resulting from applied boundary stress: tension as enlarges the wound and, in the 
elliptical case, decreases the AR. Inset illustrates relative width of thin ring, (c) Finite-element model 
detail, illustrating the mesh in the Bo configuration for the elliptical wound (circular wound is similar). 
Roller boundary conditions on two internal edges enforce quarter-symmetry conditions, and wound edges 
are stress-free. Inset details the wound region. 
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5.1. Kinematic Equations 

Applied equibiaxial tension as initially deforms the stress- free membrane Bq into the reference 



configuration bo (Fig. 5b). With erg held constant, contraction is then specified as a function of 



position and time for t > 0. The theory behind the present work is described below; further details 



can be found in previous reports on the mechanics of growth and morphogenesis (Rodriguez et al. 



1994 Taber, 2001). 




m 



/ Aa 



6(r 2 ) 



Figure 6: Sequence of configurations which transform unloaded body B into current configuration b(t), with 
fibers created at some intermediate time r. The mapping of the total deformation between times Tj and 
Tj is given by F Tj ./ r ., while G Tj / T . and F* ... give the corresponding contraction and elastic deformation, 
respectively. The contraction and elastic deformation of cells and each of the fibers is considered individually 
but all components, once created, undergo the same total deformation. 

Following the creation of a wound, cells in the thick ring begin to contract while actin and 



myosin are recruited to the wound periphery (Sonnemann and Bement, 2011 ), where they assemble 



into a thin supracellular cable that contracts over the course of minutes (Fig. [4b. As individual 
actomyosin fibers incorporate into the cable over time, they each experience a different mechanical 
environment due to the ongoing deformation of the tissue. To model this behavior, each fiber is 
assumed to form at an intermediate time r and subsequently contract. 

In the scheme illustrated in Figure [6j the first group of fibers is created at time t\ . Before 
this time, each cell in Bq contracts to a new zero-stress state defined by the contraction (growth) 
tensor G£ / , which generally varies with time and space. The cell then deforms through the 
elastic deformation gradient tensor F^*^ . This deformation is caused by surface loads, as well 
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as by geometric compatibility requirements when the cells are reassembled into the intermediate 
configuration b(r%). Fibers are created in 6(ri) at a prestretch Fq* relative to their own zero-stress 
state. Here, we assume that Fq is the same for all fibers. 

Thereafter these newly created fibers are constrained to undergo the same total deformation 
pointwise as the cells. For the time scale considered here, we assume that fiber degradation is 
negligible. During the interval from n until the next group of fibers form at T2, the cells contract 
by G^ ri , while the fibers contract by , relative to their individual zero-stress configurations, 
with the label f\ indicating the fiber under consideration. The zero-stress configurations at t± are 
obtained by reversing the elastic deformations in b(ri) of both cells (F^ ) and fibers (Fq*). The 

cells and fibers then undergo elastic deformations F^*y ri and F^t , respectively, to give the body 
6(t2) in which the second set of fibers form. This scheme continues until the current configuration 
b(t). For a constrained mixture, the deformation of the fibers, once formed, is identical to the 
deformation of the cells and the composite body as a whole, i.e., F t , = F£y r = F t / T , where we 

write f{, = F{) as the total deformation at t for fibers created at time r. 

t/T t/Ti 

The total deformation gradient tensor F t / Q maps the undeformed configuration Bq into the 
current configuration b{t). For cells, the sequence of transformations which together compose 
F t/o = F t/o is S iven W (Figure [6]) 

* t/0 - * t/r 2 • • I* r 2 /ri ) " * r 2 /n ' ^ra/n ' I* n /oJ ' * n/0 ' b T,/0 

T7C* f~iC pC pC 

= F^/ T2 • G t c /0 , (2) 
where G^ = G£, • G£ 2 , • G^ Q . This equation gives the elastic cell deformation as 

F^=F^ /0 .(G t c /0 )- 1 = Ff, (3) 

where we write F£* r2 = Ff since the cell elastic deformation is independent of the fiber creation 
time r. In general, we may write for cells {i = c) and fibers (i = f) F^ Q = F*y T • F^ Q and 
GJ/q = G|. • G^q. This decomposition does not hold, however, for the elastic deformation, i.e., 

b t/o r * t /T • *t/o- 

el 



For fibers created at n, the total deformation F^< = F t/ / T1 is (see Fig. 

t/n 



Ff 



p/l* 


t/T 2 


( F /i* ) 


T2/T1 T2/T1 V U / 


p/l* 


• g/; 

t/T2 


G /J / • 

T2/T1 


(FjT 1 


**/r 2 


t/n 







(4) 



with G w = G w • G*^ 1 / . Solving; for the elastic deformation of fibers created at n gives 



<2= F f/^ F o*-(Gf; Ti r\ (5) 
and a similar expression can be found for fibers created at t<i- In general, therefore, the elastic 
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deformation for fibers created at time r is given as 



(6) 



In the current configuration b(t), the volume fractions of the cells and fibers are 4> c and </>^, 
respectively, with (f> c + $ = 1. The fiber volume fraction in the current configuration is given by 



¥{t) J T /t dr , 



where 4>^(t) is the rate of fiber formation at time r and J T / t = JrJt 1 = det[F^ Q 
accounts for the change in fiber volume from 6(r) to b(t). 

5.2. Stress and Equilibrium 



(7) 



The total Cauchy stress tensor is given as (Humphrey and Rajagopal, 2002) 



a(t) = ^{t)a c (Ff ) + j % 0/(r) <r> (f£) J r/t dr, 



(8) 



where Ff and Ff, are given by Eqs. @ and @. We assume that the cells and fibers are nearly 
incompressible and only change shape as they contract, so that J T u = 1. 

With inertial effects being negligible, morphogenesis can be treated as quasi-static, and the 



equilibrium equation is (Humphrey and Rajagopal, 2002) 



V • a = 0, 

where V is the gradient operator defined in b(t). 



(9) 



5.3. Constitutive Relations 

Because of their high water content, soft biological tissues often are treated as incompressible 
materials. Not only does this assumption often lead to numerical challenges, but it really is not 
accurate. During deformation, water can enter or leave the tissue, as well as shift from one location 
to another within the tissue. For these reasons, we assume here that the blastoderm is nearly 



incompressible. In addition, to a first approximation, we assume that it is pseudo-elastic (Fung 



1993). 



For a compressible pseudo-elastic material, the constitutive relation can be written in the 



form (Taber, 2004) 



(10) 



i = J_ F « . 9VV 1 . ( F i*)T 

ji* dC i* y i 

Here, W l {C 1 *) is the strain-energy density function for cells or fibers, C l * = (F l *) T -F 1 * is the right 
Cauchy-Green deformation tensor relative to the current zero-stress state, and J 1 * = det F" = 1 



is the volume ratio (Taber, 2004). In the following, we drop subscripts on F and related quantities 



with the understanding that F c * depends on t while F-^* depends on t and r. 

Available experimental data suggest that tissues in the early embryo are relatively linear and 



isotropic (Xu et al. 2010; Zamir and Taber, 2004b), especially compared to mature tissues which 
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have a more organized microstructure. Hence, the cells in the blastoderm are assumed to comprise 
a compressible isotropic material whose strain-energy density function is 



W c 



2 



/l *_ 3 + l_^( (/3 * ) ^_ 1 
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which is the form for a Blatz-Ko material (Taber, 2004), where the invariants are defined as 
/* = trC c * and 1% = (J c *) 2 . In addition, [x c and v are the shear modulus and Poisson's ratio, 
respectively, in the limit of small strain. 

The fibers in the thin ring are modeled as a distinct transversely isotropic material. By con- 
struction, such fibers lie parallel to the wound margin, and the elastic stretch ratio of the fibers is 
given by 

A^* = Ve e ■ Cf* ■ eg, (12) 



incompressible neo-Hookean bar with (Taber, 2004) 



where eg is a unit vector parallel to the wound edge (see Appendix A). Each fiber is taken as an 

(13) 



A* 



+ 



2 (V* 



-i 



The fiber stress is given by Eq. ( 10 ) with 



F f* _ \f* eg e e + A T *(e r e r + e z e z 



(14) 



where the transverse fiber elastic stretch ratio X T * = 1/V \f* satisfies the incompressibility condi- 
tion detF f * = 1. 



5.4- Single-Fiber Approximation 

Equation Q is a hereditary integral which depends on the histories of deformation and fiber 
deposition. As such, it introduces significant computational cost and complexity. To simplify the 
calculation, we take advantage of the separation in time scales of the healing mechanisms. The 
first phase of wound healing, driven by cellular contraction, is essentially complete by the time 
fibers form in significant quantity, while actin staining suggests that the bulk of the thin ring forms 
within the first few minutes after wounding (Fig. [4]). Hence, cell contraction is essentially constant 
as fibers form at a rate which is relatively fast compared to subsequent contraction and healing. 
As a first approximation, therefore, we take the reference configuration for all fibers in the thin 
ring to be the same, and set r = 0. Hence, we assume that the fibers accumulate relatively rapidly 
over the first few minutes of the healing process, while subsequent contraction takes place over a 
longer time scale. Then, Eq. ^ becomes 

< = Ko = F/ *W = F f/o • H* • (Gf/o) - 1 • (15) 
With F-f* no longer dependent on r, Eq. ^ reduces to 

a(t) = (j) c (t)(T c (F c *) + a f (V*) ^ 4> f (T)dT 

= 4> C CT C + (j)f(T f , (16) 
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in which subscripts denoting time intervals have been dropped. In this approximation, the total 
deformation gradient tensor for a constrained mixture is given by F = F c = F / , and Eqs. Q and 



(15) yield 



F c * • G c = F-''* • 



(17) 



This equation gives the elastic deformation gradient tensors F c * = F-(G C ) 1 and F-^* = F • Fg* • (G-^ 
which are used in Eq. (10). With G c (t), G^(t), and Fq* specified, we employ the single fiber ap- 



proximation in the model described in the remainder of this section. In Appendix B we show that 
this approximation is sufficiently accurate for our purposes. 

5.5. Contraction Dynamics 

The cell contraction tensor G c is taken in the form 



G> e r e r + Gge e e s + G c z e z e, 



(18) 



where e r and eg are unit vectors perpendicular and parallel to the wound edge, respectively, and 
e z is normal to the plane of the membrane (see Appendix A). For convenience, we refer to the r 
and directions as radial and circumferential, respectively, even for elliptical wounds. Similarly, 
we refer to elliptical semiaxes as major and minor radii. Since a cell contracts with little change in 
volume, we take det G c = G c r G c gG c z = 1. Here, we consider three contraction schemes which define 
G%, Gg, and G c z in terms of a single specified cell contraction "stretch ratio" G c : 



Gc /~ic ric 
r , Isrg, Lr 2 



(G C ,G C , 1 

1, Lx , q, 



jQcji isotropic contraction (r, 8 directions) 



G c ,l 



circumferential contraction (9 direction) 
radial contraction (r direction). 



(19) 



As discussed below, G c is a specified function of time. Other choices of (Gy, Gg, G c z ) for circumfer- 
ential and radial contraction schemes that preserve det G = 1 , including those where G c z = 1 , do 
not significantly alter the behavior of the model (not shown). 

Since G c is isochoric and the cells nearly incompressible, contraction causes cells to thicken in 
the z direction, affecting cell stiffness (which is proportional to thickness). Because the epithelial 
membrane is restricted to lie in the plane, such thickening is the only deformation in the z direction, 
and no out-of-plane bending is permitted. 

We hypothesize that rapid cell contraction in the thick ring drives the first phase of wound 
closure. To this end, the degree of cell contraction G c varies with both position and time. Outside 
of the thick ring, we set G c = 1 at all times, indicating passive cells. Within the ring, G c decreases 
with time from an initial value of 1 to the final value G\. To avoid discontinuities, which can lead 
to numerical instabilities, the behavior of G c transitions gradually at the ring boundary. We take 
G c in the form 

G c (t,p) = l-[l-g c (t)]f(p), (20a) 
where f° and g c , respectively, govern the spatial and temporal aspects of G c : 



f C (p) 
9 C (t) 



1 I ( P~ Pc 



-i 



(l-C7f)exp(-i/T c ) + C7f. 



(20b) 
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Here, r c is the characteristic cell contraction time (within the thick ring), p c is the width of the 
thick ring, and k c determines the width of the transition at the ring boundary. The distance from 



Appendix A 


by Eqs. ( 


A.3 


) and ( 


A.6 



respectively. Figure 7a shows G c (p) for various times. 
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Figure 7: Cell and fiber contraction patterns, as well as liber volume fraction, in the vicinity of the wound 
in the undeformed configuration. Parameters are given in Tables [l] and [2} [g~ 



Cell contraction G c (Eq. 20a I 



as a function of distance from wound edge p for various times, with p (see Appendix A) normalized by 
initial circular wound diameter r,„ . Shaded region indicates domain of thick cell ring, (b) Fiber contraction 



(Eq. 22 ) as a function of time, normalized by ry . At the end of the simulation (600 sec, indicated by 
shaded region) G? has decreased from 1 to 0.46. (c) Fiber volume fraction <j>^ near the wound for various 
times (Eq. 23). Shaded region indicates nominal domain of fiber ring. Note that is spatially uniform, 
with the extent of fiber ring given by </>* . 



Consistent with the single fiber approximation, the thin contractile ring forms with a prestretch 



Fq and subsequently contracts (Gf decreases) while additional assembly occurs (<pf increases). 
Fibers are assumed to contract only along their lengths, and the fiber prestretch and contraction 
tensors, respectively, are taken as 



F f Q * = e r e r + Ffegeg + l/^F^*e z e 2 

G f = e r e r + G f e e e + l/V~Gf~e z e z , 



(21a) 
(21b) 



which satisfies det Fq 



det G* = 1, and both are symmetric relative to the cross-fiber direction. 
To define the dynamics of the fiber contraction ratio G* (t) , we assume that G* decreases expo- 
nentially with time from G$ = 1 at t = to a final value G{ at t = oo, with the rate of decrease 



given by the characteristic time constant Ty (Fig. 7b): 

Gf(t) = (l-G{)exp(-t/ Tf ) + G f v 



(22) 



The spatial distribution of the thin ring is defined through $ , which varies both spatially and 
temporally. The value of $ is uniformly zero at t = 0, increases with time in the vicinity of the 
wound at a rate characterized by r^, and approaches asymptotically the maximum value at the 
wound edge. After wounding (t > 0), the spatial distribution of <jy takes the form of a decaying 
exponential away from the wound edge whose characteristic width pj corresponds to the nominal 
thin ring width. Mathematically, 



<j) f (t,p) =<j)([l- exp(-t/"7>)] exp(-pfpf) 



(23) 
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which is plotted in Fig. 



7c 



5.6. Solution Procedure 

All model simulations were performed using COMSOL Multiphysics (v 3.5a; Comsol, Inc.), 
with plane stress conditions assumed; see Taber ( 2008 ) for more details about the implementation 
of growth. The simulations take place in two stages: starting with the undeformed configuration 
Bq, the boundary stress is increased from zero to the prescribed value <tb, where it is subsequently 
held fixed for the remainder of the simulation. There is no contraction or fiber assembly during 
this preparatory stage (i.e., G i = 1 and (j)f = everywhere). For convenience in comparing to 
experiment, where the wound is created at t = 0, we take the preparatory stage to begin at 
t = 0~ in the undeformed configuration Bq, and conclude at t = with the system in the reference 
configuration bo- The simulation stage begins at t = 0, where G c , G* , and $ vary with time as 
previously described. All results presented here are for t > 0. 

Because of the symmetry of the system, we simulate one quadrant of the epithelial membrane. 
We approximate an incision wound by an ellipse with a large AR, and use a triangular mesh of 
third-order Lagrange elements with a progressively finer mesh in the thick and thin ring regions, 



resulting in 4474 and 4419 elements in the circular and elliptical models, respectively. Figure 5c 
illustrates the mesh geometry and boundary conditions of the finite-element model. 



5.7. Parameter Values 

The model parameters can be divided into three groups. Parameters defining the geometry of 
representative circular and elliptical wounds were obtained from our experimental data. A second 
set of parameters were obtained from experiments or deduced, and are considered known. A final 
set of parameters was determined by iteratively fitting the circular wound model to experimental 
results. These parameters were then used also for the elliptical wound model. In this section, we 
discuss the choice of the first two types of parameters, which are then considered "fixed." The 
selection of the third type, considered "free" parameters, is discussed later. 



5.7.1. Geometric Parameters 

In constructing the model geometry we distinguish between the stress-free undeformed configu- 
ration Bq and the reference configuration bo, which is deformed by the boundary stress o~b but not 



fiber or cell contraction (Fig. 5b). The geometry of the model is defined in the Bo configuration, 



but is chosen such that dimensions in bo match experimental measurements. 



Geometric Parameters 


B 


bo 


Circular wound radius 
Elliptical wound major radius 
Elliptical wound minor radius 
Thick (cell) ring width 
Thin (fiber) ring width 
Membrane dimension 


38 /im (R w ) 
82 am (R a ) 

5.4 pm (R b ) 
41 pm (p c ) 

1.5 /mi (p f ) 
1020 /im 


54 pm (r w ) 
91 /im 
34 /im 
40 /im 
1.3 pm 
1140 pm 



Table 1: Wound geometries for the model in undeformed (Bo) and reference (bo) configurations (see Fig. [5]). 
Membrane dimension is the length of the square epithelial membrane quadrant (Fig. 5c I . 



The circular wound radius r w in bo was obtained from measurements immediately after wound- 



ing (see Section 4.1 ). The elliptical wound dimensions were chosen such that in bo the elliptical and 
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circular wound areas are approximately equal, and the AR (=2.67) is consistent with experiment 
(Section 4.2). Thick and thin ring widths (p c , pf) were estimated from fluorescence microscopy 
images of circular wounds (Fig. [4]), and k c is taken as /9 C /10. The same parameters are used for 
elliptical wounds. To eliminate far-field boundary effects, the dimensions of the membrane are 
chosen to be at least an order of magnitude greater than the wound. Table [T] summarizes the 



geometric parameters, and Appendix A provides details of ring geometry specification. 



5.7.2. Other Fixed Parameters 

As discussed in |Appendix C we estimate the cell shear modulus to be p c ~ 40 Pa. Literature 
suggests a wide range of fiber stiffness values, from ~ 10 kPa ( |Lu et al" 2008) to well over 1 
MPa (Deguchi et al. , 2006), with fibers much stiffer than cells (pf 3> p c ) (Rauzi and Lenne, 2011). 
We find that model sensitivity to pj decreases for pf > 100 p c as fiber dynamics alone dominate 
wound closure (data not shown), and thus choose Pf = 80p c = 3.2 kPa. We further chose 4>f to 



have a maximum value of 0.5, and note that model behavior depends on the product 
not shown.) 



fPf (analysis 



Material Parameters 


Fixed 


Cell shear modulus 


Pc 


40 Pa 


Boundary stress 




42 Pa 


Poisson's ratio 


V 


0.45 


Final fiber fraction 


4 


0.5 


Fiber shear modulus 




3200 Pa 


Free 


Cell contraction time 


T c 


12 sec 


Fiber contraction time 


T f 


600 sec 


Fiber formation time 


T 4> 


180 sec 


Cell final G 


G\ 


0.55 


Fiber prestretch 


Ft 


1.65 


Fiber final G 


G{ 


0.15 



Table 2: Cell and fiber parameters. Fixed parameters have values assumed to be known and were obtained 
as discussed in Section |5.7| Free parameter values were found by an iterative procedure so that circular 



wound model area trends reproduce those obtained by experiment (Section 6.2) 



The boundary tensile stress ob was estimated from previous work (Varner and Taber, 2010 



Varner et al. 2010). Immediately after wounding, circular epithelial wounds open to a diameter 



approximately 1.2 times that of the punch used to make them (Varner et al. 2010), a result that 
corresponds to an equibiaxial membrane stretch ratio of 1.1 (Varner and Taber 2010). We find 
that a boundary stress ob = 1-05 x p c = 42 Pa reproduces this strain. Model results are relatively 
insensitive to Poisson's ratio in the range 0.4 < v < 0.49, and we take v = 0.45. Note that nearly 
incompressible models for anisotropic materials can be susceptible to errors for v — > 0.5 (Ni Annaidh 



et al. 2012). All parameters used in the cell and fiber model, including those obtained from a fit 



to experiment (see below) are listed in Table [2j 



6. Model Results 



We now use our model and experimental data to test our main hypothesis. This requires 
determining the remaining free parameters and examining whether the model captures the main 
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features of the measured wound healing response. 



6.1. Cellular Contraction is Relatively Isotropic in Thick Ring 

Fluorescence microscopy images (Section 4.3) suggest that contraction at cell borders occurs 
within a thick ring around the wound, but do not indicate its nature. For example, is the contraction 
circumferential around the wound (like the contraction of an actomyosin cable) or isotropic? Other 
types of contraction anisotropy also may occur (Fig. 8a). Models of circular wounds are unable 
to discriminate between these possibilities, since multiple schemes can close a wound. Elliptical 
wound models, however, predict distinct temporal trends in the AR as the wound closes, and allow 
for competing hypotheses to be critically evaluated. 
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Figure 8: Area and aspect ratio (AR) trends for various cell contraction schemes for elliptical wound model 
No fibers are present in these simulations (fiber fraction <f>' = 0) and cell contraction G c (Eq. 19 ) in the ring 
is specified directly (note, in model G c < G\ = 0.55). Wound area a is normalized by area ao at t — 0. (a) 
The epithelial membrane in the cell ring can contract in both directions (isotropic) or in one direction only 
(radial or circumferential) . |bj [c]) While the wound closes for both isotropic and circumferential schemes, 
only the isotropic case yields AR trends consistent with experiment. These trends are independent of <jb- 

Consider the area and AR of an elliptical wound for three different cell contraction schemes: 
isotropic, circumferential, and radial (see Eq. (19) and Fig. [8]). As expected, radial contraction 
causes the wound to open further, allowing that scheme to be rejected. Both isotropic and circum- 
ferential contraction close the wound, but the AR trends differ qualitatively. For the circumferential 
case, the AR increases without bound, with the wound becoming increasingly slit-like as it heals. 
This behavior is at odds with experimental observations (Fig. [3]). Isotropic contraction, by con- 
trast, yields an AR that stays relatively constant as the wound closes, behavior generally consistent 
with experiment (Fig. [3]). 

We conclude, therefore, that cellular contraction in the thick ring is relatively isotropic, with 
cells on average contracting in both the radial and circumferential directions simultaneously. This 
scheme will be used for the thick ring in all subsequent models. 

6.2. Model-Predicted Changes in Wound Geometry Agree with Experimental Results 



With the geometric and material parameters established in Section 5.7 six free model parame- 
ters remain: three that quantify the degree of contraction (Fq* , G\, G{) and three that characterize 
rates (r c , r^, tj) (see Table [2]). The goal is to find a set of biologically plausible parameter values 
that reproduce the average experimental circular wound area trends, i.e., the area given by Eq. ([!]). 
Multiple circular wound simulations were performed with systematically varying parameter val- 
ues, which were adjusted iteratively to obtain an area vs. time curve that reasonably matches 
experiment. 
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Figure 9: Comparison of experimental and model wound healing behavior, (a) Circular wound. Solid curve 
is the experimental wound area given by Eq. ([I]) with mean parameters from Table SI Circles indicate 
model results (parameters from Table [2]). (b,c) Elliptical wound. Solid curves show experimental area (b) 
and AR (c). Model results (dashed curve) use same parameters as circular wound. Insets detail model 
behavior in first minute, with same axis labels as main plot. 



The healing response during the first 20-30 seconds is dominated primarily by cell contraction, 
as the fibers have not yet formed in significant quantity. The wound behavior during this time is 



controlled by the rate of cell contraction r c and the final cell contraction parameter G\ (see Eq. 20b). 



The values G\ = 0.55 and r c = 12 sec reproduce this initial stage of contraction (Fig. 9a). 

The parameters corresponding to the fiber prestretch (Fq*) and contraction (G{) govern the 
remainder of wound closure. We find that fibers forming rapidly (r^ = 180 sec) in a state of 
significant prestretch (Fq* = 1.65) and thereafter contracting more slowly (r/ = 600 sec, G{ = 
0.15), together with the other parameters in Table [2l reproduce the experimentally observed area 



versus time curves quite well (Figure 9a ) . The uniqueness and biological plausibility of these values 
are discussed later. 

Once obtained for circular wounds, these same parameters were then used in the elliptical wound 
model, with the resulting area and AR plotted in Figures 9b, c} From an experimental perspective, 



elliptical wounds are less reproducible, and we restrict ourselves to a qualitative comparison of 
wound trends. As in the case of circular wounds, and consistent with experimental trends, the 
elliptical wounds in the model close in two phases - a rapid initial closing followed by slower 
healing (Fig. |9b[ ). The AR displays a brief decrease followed by a rebound and relatively constant 
value (Fig. |9c[). Both features are broadly consistent with our experimental observations. 

6.3. Peak Stresses Increase as Wounds Close 

Stress distributions near the wound boundary are shown for circular wounds, with radial (a r ) 



and circumferential (o~g) stress components plotted for four illustrative time points (Fig. 10). At 
t = (i.e., bo configuration), the membrane is stretched passively by o~b- Significant cell contraction 
has occurred by t = 12 sec, but fiber formation is negligible. At t = 300 and 600 sec, cell contraction 
is fully developed, with changes in stress driven by ongoing fiber formation and contraction. 



For a circular wound, a r develops two peaks as the wound closes (Fig. 10a). One peak occurs 
near the wound edge, while the other develops just outside the thick ring. The distribution of oq 
is dominated by a strong peak at the edge of the wound, increasing from about 2<tb at t = to 
40 a b at t = 300 sec. Thereafter oq decreases at the wound edge and becomes slightly compressive 
immediately outside the cell ring; the reason for this effect is discussed later. In all cases, both 
stresses approach the far-field value o~b away from the wound. 
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Figure 10: Epithelial membrane stress in the vicinity of circular wound. Icon indicates stress component 
(radial oy, circumferential erg). Distance from wound edge in the reference configuration is given in multiples 
of the circular wound radius r w , and the fiber and cell regions are indicated with blue and red shading, 
respectively (see Fig. [7]). Inset in (b) details stress near wound boundary, with axes labels as in main figure. 
Far from the wound, both stress components approach the far-field value Ob = 42 Pa. 



Stress distributions near an elliptical wound exhibit similar characteristics (Fig. S4). For an 



ellipse, the strongest stress concentrations occur near the wound edge along the major axis, while 
the stresses are reduced relative to those for a circular wound along the minor axis. 

6.4- Elliptical Wounds Briefly Open at Low Prestress 

The principal effect of varying <tb is to change the initial size of the wound. Also, as as 



increases, elliptical wounds become initially rounder (Fig. 11a). 



For a circular wound, area decreases in a similar manner regardless of o~b, with wounds under 



higher tension closing more rapidly (Fig. lib). Elliptical wounds at high tension likewise begin 
to close immediately after contraction begins, but for gb < 15 Pa the wound area increases 
momentarily before closing (Fig. lib). This initial "gaping" of elliptical wounds at low o~b is 



associated with a rapid decrease in the wound AR (Fig. 11c). Similar behavior was observed 



experimentally, where some elliptical wounds initially increased in area while becoming more round 
(Fig.§. 
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Figure 11: Effect of prestress <7b on circular and elliptical model wounds, (a) Wound size and shape at t = 
(ao) for varying as- Prestress is relative to ob = 42 Pa (Table [2]). (b) Area for circular wounds (dashed 
curves) and elliptical wounds (solid curves) with varying as- (c) Aspect ratio (AR) for elliptical wounds. 
Low initial tension in elliptical wounds results in a brief increase in wound area and a sharp decrease in AR. 



This transient area increase is driven primarily by a brief expansion of the minor radius of the 



wound (Fig. S5a). This effect is captured by the elliptical wound model when the initial tension 
is reduced. With ob set to the estimated physiological value o~b = 42 Pa (Table [2]), the area 
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decreases immediately after wounding (Fig. llb[ c) and both the major and minor axes decrease 
(Fig.[S5bl), with the AR remaining relatively constant (Fig. |1 lcj ) . For low as wound area increases, 
the minor axis increases, and the AR decreases briefly just after contraction begins (Fig. llb|c 



S5b). The reason behind this behavior is discussed later. The qualitative aspects of wound closure 



are relatively insensitive to wound size and shape (see Supplemental Appendix S.3) 



7. Discussion 

Our results suggest that wounds in the early chick embryo heal by three distinct and essentially 
sequential mechanisms. First, during about the first 30 seconds, a ring of cells (3-4 cells deep) 
contracts to quickly close the wound area by more than 50%. Second, a relatively thin contractile 
ring forms at the wound edge and contracts more slowly to close the wound nearly completely over 
a period of several minutes. Finally, filopodia pull and zip the edges of the wound together to 



complete the healing process. The last two phases have been known for some time (Jacinto et al. 



2001 Woolley and Martin 2000), and the initial rapid healing phase was observed in chick embryos 



by Bortier et al. (1993), but they did not investigate it in detail. Notably, this phase seems to be 



2011 


Hutson et al. 


2009 


Ma et al. 



2009), so it may not be a universal phenomenon. As discussed below, our model captures the 
fundamental behavior of the first two phases of wound healing and yields insight into the detailed 
mechanics. 



7.1. Parameter Values for Contraction Dynamics 

For our wound healing model to be plausible, it is important that the values of the free parame- 
ters are consistent with those reported in the literature for related systems. Our model incorporates 
three separate time constants (Table [2]): r c for cell contraction, for fiber formation, and tj for 
fiber contraction. Three additional parameters specify the degree of initial and final contraction: 
final cell contraction is given by G\, and initial (prestretch) and final fiber contractions are given 
by (Fq*) _1 and G{, respectively. 

The rate of cell contraction r c (12 sec) matches the time constant t\ (12.4 sec) for the first 



phase of healing circular wounds (Section 4.1 Table SI). Both the speed and magnitude of cell 



contraction (G\ = 0.55) are in line with observations for smooth muscle cells (An and Fredberg 



2007). The rate of fiber formation ta, (3 min) is consistent with our observed time rate of fiber 
formation (see Fig. [4] and Sec. 7.2.1), and the fiber contraction rate tj (20 min) is consistent with 



actomyosin fiber dynamics continuing over tens of minutes following a stimulus ( An and Fredberg 
2007 Brock et al.[|1996[). Embryonic t issues can shorten by 70% or more within a relatively short 
time period (Varner and Taber 2012b), supporting the value G{ = 0.15. Such a large contraction 
likely involves significant remodeling of actomyosin fibers, similar to that observed in smooth 



by Martin et al. (2009) 



muscle cells (Matsumoto and Nagayama, 2012) or the contractile ratcheting mechanism described 



Fina 



ly, the dynamics of wound closure are quite sensitive to the value of the prestretch Fq* 
(Fig S3a). Our value {Fq* = 1.65) is somewhat larger than that reported in the literature (1.10 to 
1.35; Kaunas and Deguchi 2011). This is discussed further in Section 7.2.2 

The parameters in Tableware not unique. For instance, a set of parameters where fibers form 
slowly (tj = 20 min) with a larger prestretch Fq* = 3.1 and without further contraction (G{ = 1) 



can also reproduce our wound closing results (Fig. S3b). Our experimental data are insufficient to 
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discriminate between these and other possible parameter sets, and this certainly warrants further 
study. 



7. 2. Phases of Embryonic Wound Healing 



Unlike adult wounds, embryonic wounds generally heal quickly without leaving a scar (Redd 



et al. 2004). The healing process involves many of the same mechanisms used for morphogene- 

;issue shape driven 
during neurulation 



sis (Wood et al. 2002), which often produce dramatic and rapid changes in tissue shape driven 
by contractile forces (Davies, 2005). Fusion of epithelia is also common, e.g. 



and heart tube formation (Colas and Schoenwolf, 2001 Moreno-Rodriguez et al. 2006, Ray and 



Niswander, 2012). The three healing phases close wounds efficiently and robustly. 



7.2.1. Phase 1: Rapid Cellular Contraction (Thick Ring) 

Activation of upstream GTPase regulators of actin and myosin occurs 10-20 seconds after 



wounding (Bement et al. 2006 Clark et al. 2009), and wound-induced actomyosin structures are 



first observed 1-3 minutes following wounding in a variety of organisms (Brock et al. 1996; Kasza 



and Zallen, 2011). These characteristics are consistent with our observations on the formation and 



contraction of the thin contractile ring (see Sec. 7.2.2 and Figs. 4b|e ). However, we have found that 
healing actually begins within seconds of wounding, ruling out the possibility of de novo actomyosin 
assembly in the first phase. 

Presumably, the initial rapid-healing phase is instead powered by pre-existing cellular cytoskele- 
tal structures. Epithelial cells generally contain a ring of actin and myosin associated with adherens 
junctions around cell borders (Lodish et al. 2004) (Figs. 4a|d ), and contraction of these fibers 
plays a role in many morphogenetic processes (Martin, 2010 Rauzi and Lenne, 2011). Shortening 



of these relatively randomly oriented fibers would result in tissue-level isotropic contraction, con- 
sistent with our model. A possible signal for the contractile response is calcium (Ca 2+ ), which can 



induce swift but short-lived actomyosin contraction (Cordeiro and Jacinto 2013, Tomasek et al. 



2002), is rapidly modulated (Clark et al. 



cell layers away from the wound (Klepeis et al. 



2009 



Xu and Chisholm, 2011), can be activated several 



2001 Woolley and Martin, 2000 Xu and Chisholm 



2011 ), and has been implicated in a variety of wound healing processes (Benink and Bement , 2005 



Sonnemann and Bement, 2011). Thus, the proposed cell contraction mechanism is reasonable from 



a cellular biological perspective. 

Most previous studies of chick embryonic wound healing have focused on time scales of a minute 
or longer, but Bortier et al. (1993) have described rapid dynamics of chick epithelial wounds 



similar to those we report here. They did not, however, comment on its significance. The authors 
note that "the wounds showed a viscoelastic reaction upon wounding: they enlarged within one 
second, immediately followed by a narrowing, until wound diameters were between the original and 
largest size. After these reactions the wound submarginal region was thickened..." It seems likely 
that the submarginal region, reported as 75-150 fim wide, corresponds to the thick ring where 
isotropic contraction takes place. The authors also observed blebbing cells in this region, while 
cell bulging and "mounding", as well as tissue thickening, has been reported near the margin of 



embryonic wounds in several organisms (England and Cowper, 1977; Jacinto et al. 2001; Smedley 
and Stanisstreet 1984, Stanisstreet et al.] ~1980). Both cell bulging and blebbing suggest high 
intracellular pressure ( Charras et al. 2008 ) , which is consistent with strong isotropic contraction. 
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7.2.2. Phase 2: Slow Fiber Contraction (Thin Ring) 

Following the creation of a multicellular wound, actin and myosin are recruited to the wound 
margin where they assemble a continuous supracellular cable which spans cells through adherens 
junctions. The cable forms over the course of minutes and contracts to generate tension which 



helps close the wound (Brock et al. 1996 Clark et al. 2009 Kiehart, 1999 Martin and Lewis 



1992 Sonnemann and Bement 2011). Contractile actomyosin cables are used by the embryo in 



a variety of contexts, for example to divide cells and close the dorsal ectoderm in Drosophila, as 



well as to heal wounds (Davies, 2005 Rodriguez-Diaz et al. 2008 Sonnemann and Bement, 2011). 



Studies of actomyosin dynamics at high temporal and spatial resolution have shown contraction 
to be pulsatile, with periods of active shortening occurring between periods of quiescent stabi- 
lization (Gorfinkiel and Blanchard 2011 Levayer and Lecuit 2012; Martin 2010). We did not 



observe such pulsatility, as localized pulses may be averaged out over the perimeter of a multi- 
cellular wound, and we did not include this in the model. Being relatively stiff, the fiber ring can 
generate the force needed to overcome the high circumferential stresses that develop as the wound 



closes (Fig. 10) (Strieker et al. 2010). 



The remodeling theory used in our model was originally developed to describe the turnover 
of passive tissue constituents, e.g., elastin and collagen, which are thought to be created with 



significant prestretch (Humphrey and Rajagopal, 2002). It is not clear, however, that actomyosin 



fiber assembly involves an initial stretch. Hence, we interpret Fq* as the result of a rapid contraction 
of the initially unloaded fiber (time constant r un <C 7/) that slows (r un — > t/) as tension increases, 
with the observed prestretch Fq* = (G^) _1 resulting from the contribution of many such fibers. 



This Hill-like behavior is consistent with some recent models for stress fibers (Deshpande et al. 



2007 Stachowiak and O'Shaughnessy 2008), and an ability of actomyosin fibers to rapidly contract 



by large amounts has been observed in smooth muscle cells (An and Fredberg 2007). 



7.2.3. Phase 3: Filopodial Zippering 

Filopodia-mediated zippering is also an important wound healing mechanism in early chick 
embryos. Although cell protrusions have not been reported during wound healing in older chick 



been seen at early (HH 3-5) embryonic stages (England and Cowper 



embryos (Brock et al. 1996| Martin and Lewis, 1992), filopodia, lamellipodia, and microvilli have 



Stanisstreet et al. , 1980[ ), and were apparent here as well (Fig. |4f SI ). Such protrusions can contact 



1977 Mareel and Vakaet , 1977 



and pull cells at opposite sides of the wound together, closing it by a zippering mechanism (Jacinto 
Wood et alT| |2002[ ). Because of their short length (~5 /im ( |Wood et alj |2002[ )) 



et al., 2000 



filopodia and related structures are unlikely to play a significant role until the wound opening is 
relatively small. Such effects are not captured by our model. 

The three active processes - isotropic cell contraction, actomyosin ring contraction, and filopo- 
dial zippering - work together to close a wound, a redundancy which has been observed in other 
morphogenetic processes (Davidson et al. 20021 iJacinto et al. 2001 Kiehart et al. 2000). Very 



small wounds tend to close quickly (data not shown), with isotropic contraction likely playing the 
dominant role. The capacity of this mechanism is limited, however, and large wounds require a 
contractile ring. In some cases wounds never heal, or heal so slowly that closing does not occur 
during the period of observation. The reason may be that they are created in areas of older embryos 



where blastoderm tension is high (Varner et al. 2010), or they excessively perturb the tissue past 



the point of viability. Nevertheless, the capacity of the early chick embryo to heal is remarkable, 
and we find that nearly all wounds, if given enough time, do close. 
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7.3. Mechanics of Wound Closure 

The way circumferential contraction closes a wound is consistent with physical intuition, i.e., 
by shortening the perimeter. The mechanics of embryonic wound healing, however, contain some 
subtleties that warrant discussion. 

For example, it is important to note that cellular contraction occurs with little or no change in 
cell volume. Hence, circumferential contraction is accompanied by radial expansion, which helps 
healing by pushing the edges of the wound inward. Because of structural stiffening in regions of 
high curvature, these radial stresses deform edges along the minor axis of elliptical wounds more 
easily than those along the major axis, causing the wound to become more slit-like. This explains 
in part why the AR increases for the case of circumferential contraction only (Fig. 8c). In contrast, 
radial contraction pulls the long edges outward, increasing the AR (Fig. 8c). 

Since radial and circumferential contraction in the thick ring have opposite effects on the AR, 
these two effects essentially cancel out for isotropic contraction, and the wound shape remains 
relatively constant as it closes (Fig. 8c). The interaction of these effects with stress-induced elastic 
deformation underlies the initial opening of the elliptical wounds in membranes subjected to rela- 
tively small tension (Fig. lib). As the applied tension increases, the stress concentrations near the 
edge of the wound increase. Since radial stress is considerably smaller than circumferential stress 
near the wound edge (Fig. S4), radial contraction is more effective initially than circumferential 
contraction, and the wound first opens before starting to close. Hence, the detailed wound geom- 
etry depends on a balance between radial and circumferential contraction, as well as the tissue 
stresses developed in these directions. 

Another item of note is that inward radial movement of the membrane reduces its perimeter 
and tends to decrease circumferential stress. Unless this "slack" is taken up by circumferential 
contraction, compression can result. This effect is clearly seen in plots of o~g (Fig. 10), where the 
circumferential stress decreases, and in some cases becomes negative, immediately outside of the 
thick ring, as well as inside the ring toward the end of the simulation. In a physical system such 
compression, if large enough, would manifest as buckling or wrinkling of tissue, which we have not 
observed experimentally. 

The stress concentrations predicted by our model, particularly around elliptical wounds, are 



of a magnitude consistent with fiber stress exerted by smooth muscle cells (An and Fredberg 



2007). Viscoelastic effects, not incorporated in the model, are likely to reduce the peak stress 
considerably. Overall, the coupling between passive stretch, active contraction, and nonlinear 
geometric effects driven by large tissue deformation resist simple explanation and highlight the 
utility of computational models. 

Mechanical signals govern the spatiotemporal contractility of tissue in a variety of ways ( Gorfinki il 



and Blanchard] [20TT| |Hutson and"l4a] [2008} |Levayer and Lecuit| [20T2| |Vogel and Sheetzj [2006} 



Wozniak and Chen 



2009), and mechanical feedback has long been suggested to play a role in 



wound healing (Martin and Lewis, 1992). While we have recently proposed a mechanical feedback 
model for embryonic wound healing (Taber, 2009), this is not included in the current model. 

Finally, it is important to note that other mechanisms may be involved in embryonic wound 
healing. Some of these may be redundant mechanisms that activate when contraction fails. For 
example, increased adhesion affinity between cells can cause them to elongate radially and shorten 
circumferentially, creating tension that helps close the wound. Hence, several aspects of this 
problem warrant additional future study. 
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Appendix A. Contractile Ring Geometry 



For clarity, in this appendix we use the notation (r, 9, z) for cylindrical polar and (£, r/, z) for 
elliptical cylindrical coordinates, respectively. In the body of the manuscript, for convenience we 
use (r, 6, z) to refer to coordinates in both coordinate systems, with the understanding that £ = r 
and 77 = 9 in the context of an elliptical geometry. 

Appendix A.l. Coordinate Systems 

Because the contractile fibers are oriented parallel to the wound edge, we construct coordinate 
systems with a unit vector lying parallel to the fibers in the undeformed state Bq. For circular 
wounds we use a cylindrical polar coordinate system, which relates the coordinates (r, 9, z) to the 



Cartesian (x,y,z) coordinates as (Fig. A. 12a) 



x = rcos9, y = rsin9, z = z (cylindrical). (A.l) 
For an elliptical wound centered at the origin with foci at positions ±a on the x axis, the coor- 



dinates (x,y,z) relate to the elliptical cylinder coordinates (£,77,2) as (Arfken and Weber, 2005) 



(Fig. A.12b) 



x = a cosh £ cos??, y = asinh£sinr/, z = z (elliptical). (A. 2) 

The unit vectors e r and e^ are normal to the wound edge in the polar and elliptical coordinate 
systems, respectively, while eg and e^ are perpendicular to the wound edge. In both cases e z is 
normal to the plane of the membrane. 




Figure A. 12: Detail of model geometry in the immediate vicinity of the wound in the undeformed (-Bo) 
configuration for circular (a) and elliptical |(b)| wounds. Because of symmetry, only one quadrant of the 
membrane is considered, and the wound lies in the lower-left corner. The solid red curve indicates the 
wound boundaries, and the green curves mark the boundaries of the thin and thick rings. Dashed blue 
and red curves indicate direction along the circumferential (eg and e^) and radial (e r and e^) directions, 
respectively. The vector p marks the distance between an arbitrary point X and a corresponding point R 
on the wound edge. 



Appendix A. 2. Contractile Ring Width 

To construct rings of uniform width around a wound in the undeformed state Bq , we calculate 
the distance from any given point X to the wound edge. For a circular wound this task is straight- 
forward. We define R as the nearest point to X on the wound edge, and compute the magnitude 
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of the vector p = X — R as (Fig. A. 12a) 

p = rx — rR (cylindrical). (A.3) 
An analogous approach fails in elliptical coordinates because the quantity £x — £r does not 



correspond to a uniform distance around the perimeter of the wound (Fig. A. 12b). Instead, we 
define R as the point on the wound edge with the same r/ coordinate as X, and expand X as a 
Taylor series in = £x — £r about R. Retaining terms to 0((f£ 2 ), we write 



= X - 


R, 






d~R 








dR 




1 a 2 R 


" ~w 





1 d 2 K 



2 



R 



R, 



df. 



(A.4) 



R 



Using Eq. (A.2), we can write the derivatives of R = xex + yey a s 

= asinh^cosr/xex + acosh^sinr/xey, 



o 2 k 



R 



R 



ot cosh ^r cos fjx^x + asinh^/j smrjxey- 



(A.5) 



Finally, from Eq.(A.4) we find the magnitude of p for the elliptical wound as 

[(2 sinh £,ndl; + cosh ^nd^ 2 ) 2 cos 2 r\x + 



a 
2 



(2cosh£ft(i£ + sinh£#<i£ 2 ) 2 sin 14 r\x I (elliptical), 



1 1/2 



(A.6) 



where £r is constant along the wound edge and ijx is given by X. While this approximation strictly 
holds only for X in the immediate neighborhood of R, it serves as an adequate approximation over 
the domain of both contractile regions. 

Given the model geometry (Table [T]), we find the focus position a = 0.0798 and £ at wound 
boundary £ r = 0.0663. The contours of p equal to pj and p c , which demarcate the borders of the 



fiber and cell contractile regions, respectively, are plotted in Fig. A. 12b and we find that p in the 
elliptical coordinates evaluated along the X and Y axes in the Bq configuration varies less than 



Appendix B. Comparison of Single and Multiple Fiber Models 

To evaluate the accuracy of the single-fiber approximation, we consider the simple case of an 
unstretched homogeneous bar fixed at both ends (Fig. B.13 inset), and neglect the contribution of 
cells by setting p c = 0. In this one-dimensional approximation, F?* = A^*, Fq* = Aq*, C^* = A^* , 
and Jf* = 1. For simplicity we set p,f = 1 and Aq* = 1, and write the stress as a linear function 
of the elastic stretch ratio, 

a f = \f* - 1. (B.l) 
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For the multiple-fiber model, the fiber stretch ratio is given by Eq. @ MS 



A 



/* 

t/r 



t/T 



G 



(B.2) 



and for the single-fiber model, Eq. (17) gives 



t/r 



(B.3) 



i f 

In the present isometric problem, the total stretch ratios are X t , = A = 1. 

With the above equations, the fiber stress, given by the 1-D forms of Eqs. ^ and (16), is 

1 



OMF{t) 

for the multiple-fiber model and 

&SF 



G 



f 

t/T 



1 dr 



(B.4) 



(t) = <P f {t) 



1 

Gf 



(B.5) 



for the single-fiber model. 

For illustrative purposes, we take G$ (i) = exp(— t/r A for the single-fiber model. This gives 



G 



f 

t/T 



Gf(t)/Gf(r) = eM-(t-r)/r f ] 



for the multi- fiber model. The fiber fraction, from Eq. (23), is 

<fS(t) = l-exp(-t/r^). 
We also nondimensionalize time by the fiber contraction rate and define 

t' = t/Tf, 

a = Tf/r^. 



(B.6) 



(B.7) 



(Bi 



Finally, we can write Eqs. (B.4) and (B.5) as 



exp(— at') + exp(t') 



VMF(t') = i , . 

1 + a 1 + a 

<r SF (t') = [1 - exp(-ai')] [exp(t') - 1] 



(B.9) 



where the multi-fiber stress was integrated analytically. Both stresses are parameterized by a, 
which is large when fibers form quickly and contract slowly, and small when the opposite is true. 

The stresses omf and gsf are plotted for a range of a (Fig. B.13 ). Differences between the single 
and multi-fiber models decrease with increasing a. As a — > oo all fibers appear instantaneously 
at t = 0, so that the two models yield identical results. The parameters used in the body of the 
manuscript correspond to a = 3.3 (Table We conclude that the single fiber model, while a 
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1,2 3 

t 



Figure B.13: Fiber stress in a homogeneous bar with fixed ends (inset) as a function of non-dimensional 
time t' for different fiber models and values of a. In all cases asF > <^mf, but this difference decreases with 
large a. 



simplification of a more biologically realistic multi-fiber formulation, is sufficiently accurate for the 
chosen parameters to justify its significant computational advantages. 



Appendix C. Mechanical Properties of Blastoderm 

To estimate the mechanical properties of the early chick blastoderm (i.e., a value for fi c ), we 
constructed a simple linear model for indentation of a plate under in-plane tension, and used it to 
simulate microindentation experiments of HH stage 5 embryos previously conducted in our labora- 



tory (Varner et al. 



2010 



Varner and Taber, 2012a). Briefly, following the framework employed by 



(Zamir and Taber, 2004a), we model the blastoderm near the indenter as a thin, isotropic annu- 
lar plate under in-plane tension. For small deflection, the governing differential equation is given 



by (Szilard, 1974) 



T r V 2 w=p, (C.l) 

where w is the transverse deflection, V is the gradient operator, D p = Eh 3 /12(1 — u) 2 is the flexural 
rigidity, T r is the radial in-plane force per unit length, p is the applied surface pressure, E is the 
Young's modulus, v is Poisson's ratio, and h is the plate thickness. 

Here, we consider the case of a circular plate clamped at the outer radius b, which represents 
the distance at which the local disturbance from the indenter is effectively zero. The plate also is 
taken as fixed at the inner radius to a rigid cylindrical indenter of radius a that exerts a force P. 
The boundary conditions take the form 



V(a) = - 

9(a) = 
w(b) = 
0(b) = 



P 

2na ' 



(C.2a) 



where V is the transverse shear per unit length and O is the rotation. Note that, in this case, no 
surface pressure is applied to the plate (i.e., p = 0), so the applied load enters not through any 



specific term in equation (C.l), but rather through enforcement of the boundary condition (C.2). 



Hole punching experiments in HH stage 5 embryos have indicated that the blastoderm is initially 



in an approximately uniform, equibiaxial state of tension (Varner and Taber 2010; Varner et al. 
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2010). In addition, comparing the geometry of the wounds to the dimensions of the punching 



pipette have revealed that the initial tissue strains are e r = eg = e ~ 0.1. For plane stress, the 
prestress is then given by 



E 



O r = Og 



1-V 2 

Ee 



1 - v 

If we assume material incompressibility (i.e., v = 0.5) and e = 0.1, then 

T r = ah = 0.2Eh. 



(C.3) 



(C.4) 



Microindentation experiments performed in our laboratory have indicated an approximate bias- 
todermal stiffness of 0.2-0.6 mdyne//mi using an indenter with a radius a = 10 fim (Varner et al. 



2010 Varner and Taber, 2012a). Moreover, optical coherence tomography reconstructions of early 
chick embryos have shown the thickness of the blastoderm to be approximately h = 50-90 /im 



(Varner et al., 2010) 



If we use these values in combination with equation (C.4) and take the distance from the 
embryonic midline to the boundary of the area pellucida (~l-2 mm) as the radius of the circular 



plate b, we can solve equations (C.l) and (C.2) for different values of E and compare the computed 



stiffness values, given by P/w(a), to those reported experimentally ( Varner and Taber[ |2012a ). 
Doing so indicates an approximate Young's modulus of 120 Pa for the blastoderm, which (given 
material incompressibility) yields the shear modulus fi c = E/3 = 40 Pa. 
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Appendix S. Supplemental Material 

Appendix S.l. Parameters for Experimental Wound Area 



Co [^m 2 ] 
ci [fim 2 ] 

Ti [sec] 
t-2 [sec] 



1 



-160 
7090 
3000 
13.5 
186 



60.2 
5320 
2440 
8.83 
82.7 



459 2110 

4170 5220 

1780 3350 

9.62 16.1 

111 131 



199 
4690 
2430 
12.4 

134 



6 



230 
8850 
3550 
13.7 
195 



Mean ± SD 



483 ± 751 
5890 ± 1600 
2760 ± 606 
12.4 ± 2.5 
140 ± 40 



Table SI: Parameters for best fit (Eq.[T]) to experiment for six circular wounds. Numbers in column header 
indicate the wound in Fig. la colors match the plot in Fig. [2] for each wound. 



Table SI lists the parameters obtained from a best fit of Eq. ([!]) to closing trends of six circular 
wounds (Figs. [I] and [2J. Mean parameter values yield area versus time curve against which model 
results are compared. 

Appendix S.2. Actin and Myosin Staining 




Figure SI: Detail of actin staining in Fig. 4e with cells (gray), nuclei (blue), and actin (green). Arrows 
indicate lamellipodia; scale bar 50 urn. 



Fluorescent staining for actin at one minute after wounding (Fig. SI) illustrates the formation 
of the actin purse string and the presence of lamellipodia. Staining for phosphorylated myosin 



light chain (pMLC) at different times after wounding (Fig. S2 ) yields distributions similar to those 



for actin (Fig. [4]). These results indicate that pMLC co-localizes with actin in the thick and thin 
rings, and demonstrate that these structures are actively contractile. 

Appendix S.3. Parameter Sensitivity Studies 

The dynamics of wound closure are quite sensitive to the value of the prestretch Fq* , as illus- 
trated in Fig S3a Parameter sets other than those listed in Table [2] can reproduce the experimental 



circular wound trends, as shown in Fig. S3b 



The principal difference between circular and elliptical wounds is more intense stress concen- 
trations along the major axis of the ellipse (Fig. S4). Along both axes, at t = the radial stress a r 



is zero at the wound boundary and increases with distance, converging to o~b some distance away. 
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10 min 



(b) 



Figure S2: Fluorescence images of ectoderm near circular wound stained for phosphorylated myosin light 
chain (pMLC) at given times post wounding. Scale bar 100 \xm. 
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Figure S3: Effects of parameter variation on model behavior, (a) Effect of prestretch Fq* on circular 
wounds, with all other parameters as in Table [2J Fq* has a strong effect on wound healing behavior, 
(b) Comparison circular wound experiment and model results for two different parameter sets. "Default" 
parameters correspond to those in Table [2J For "Constant G?" parameters, G[ = 1 (no fiber contraction) 
with = 20 min and Fq* — 3.1. Results are essentially indistinguishable given the natural variation in 
wound behavior. 
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Figure S4: Epithelial membrane stress in the vicinity of elliptical wound. Icons indicate position (along 
minor or major axis for elliptical wound), and stress component (radial oy, circumferential erg). Distance 
from wound edge in the reference configuration is given in multiples of the circular wound radius r w , and 
the fiber and cell regions are indicated with blue and red shading, respectively (see Fig. [7]). Far from the 
wound both stress components approach the far- field value <jb — 42 Pa, indicated by black dashed line. 



However, the increase is very sharp along the major axis and much more gradual along the minor 
one. Similarly, the circumferential stress at t = has a value about 70 times greater than as at 
the wound border along the major axis, while it decreases to less than half of <tb along the minor 
axis. Otherwise, contraction of cells and fibers changes the stress distribution in a manner broadly 
similar to that observed for the circular wound (Sec. 6.3). 

Differences in membrane pre-stress can explain the initial variations in the healing behavior 
observed in elliptical wounds. Fig. [S5] plots the major and minor radii of both experimental and 
model elliptical wounds. The initial opening of elliptical wounds is driven by an increase in the 
minor radius, a behavior which is reproduced by the model for low values of <7g (see Sec. 6.4). 

To evaluate the role of wound size and shape on the rate of wound closure, we varied the 
initial wound size and aspect ratio. With other parameters unchanged, a smaller wound implies 
proportionately wider cell and fiber contractile rings. The model predicts that smaller wounds 
close more quickly, particularly in the initial phase of healing where cell contraction dominates 
(Fig. S6a). The aspect ratios of smaller wounds increase somewhat more than those of larger 
wounds in the second, fiber-driven phase of wound healing (Fig. S6b), an effect which is analogous 



with the result of Fig. 8c In general, however, qualitative aspects of wound closure tend to be 
relatively insensitive to wound size. This is broadly consistent with experimental observations, 
where we find that wounds much larger than those considered here will typically heal, albeit much 
more slowly. 

Changing the elliptical wound geometry by modifying the minor radius yields negligible changes 
in the relative rate of wound closure, and predictable variation in the aspect ratio (Fig. S7). Thus, 
neither the wound size nor shape has a strong effect on wound behavior, and at least for the value 
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Figure S5: Major and minor radii of elliptical wounds in initial stages of closing. Solid and dashed curves 
represent major and minor radii, respectively, as indicated by ellipse icons. |(a)| Radii of experimental 
wounds. Colors represent individual experiments (data from Fig. [3]). |(b)| Radii of model wounds, with colors 
indicating boundary tension <tb, data from Fig. 11 
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Figure S6: Wound area and aspect ratio (AR) with varying initial area (same ob)- Dashed and solid curves, 
correspond to circular and elliptical wounds, respectively, as indicated by icons. Undeformed wound radii 
(R w for circular wound, Ra and Rg for elliptical) are varied relative to nominal values (Table[TJ). Decreasing 
wound size increases effective width of fiber and cell rings, and results in slightly faster wound closing. 
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Figure S7: Area and aspect ratio (AR) of model elliptical wounds with varying initial AR. The minor radius 
in undeformed configuration R% is varied relative to nominal value (Table [l|, i.e., Rb — 0.5Rb specifies a 
wound with increased AR, and a circular wound is equivalent to Rb — 15Rb- Wound geometry variation is 
reflected mainly in initial AR and has little effect on relative area or AR trends. 
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of gb under consideration, circular and elliptical wounds behave similarly. 
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